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Abstract 

We show how Newtonian gravity emerges on 4-dimensional non-commutative 
spacetime branes in Yang-Mills matrix models. Large matter clusters such 
as galaxies are embedded in large-scale harmonic deformations of the space- 
time brane, which screen gravity for long distances. On shorter scales, the local 
matter distribution reproduces Newtonian gravity via local deformations of the 
brane and its metric. The harmonic "gravity bag" acts as a halo with effective 
positive energy density. This leads in particular to a significant enhancement 
of the orbital velocities around galaxies at large distances compared with the 
Newtonian case, before dropping to zero as the geometry merges with a Milne- 
hke cosmology. Besides these "harmonic" solutions, there is another class of 
solutions which is more similar to Einstein gravity. Thus the IKKT model 
provides an accessible candidate for a quantum theory of gravity. 
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1 Introduction 

The aim of this paper is to study the physical properties of the effective ("emergent") 
gravity which arises in matrix models such as the IKKT model, with emphasis on long- 
distance and cosmological aspects. 

The presently accepted description for cosmology, known as ACDM model, is able to 
reproduce the basic observational data. Its is based on the assumption that gravity is 
described by general relativity (GR), leading to the Friedmann equations which describe 
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a homogeneous and isotropic universe. However, this requires a very dehcate fine-tuning 
of the basic parameters in order to reproduce the observational constraints. The most 
problematic aspect is dark energy or the cosmological constant A, which must be fine- 
tuned to A ^ 2 X lO^^eV on order to produce the apparent cosmic acceleration inferred 
from the type la supernovae data. The standard model also requires significant amounts 
of dark matter in order to reconcile the galactic rotation curves with GR: The rotation 
velocities of stars near or outside the visible bulk of typical galaxies does not decrease with 
radius as implied by GR resp. Newtonian gravity, rather it is more-or-less fiat or even 
increasing. This is usually explained by postulating large dark matter halos around the 
galaxies, which however has never been observed directly. 

These and other problems provide sufficient motivation to look for alternative descrip- 
tions of gravity, in particular at large distances. Indeed in spite of its undisputed success, 
general relativity has never been tested directly on cosmological scales. Finally, all at- 
tempts to define a quantum version of GR are faced with serious conceptual and technical 
difficulties. 

In this paper, we show that matrix models, notably the IKKT matrix model [1], provide 
a very interesting alternative theory of gravity which might allow to resolve these problems. 
The geometric mechanism for gravity in these Yang-Mills matrix models was clarified re- 
cently [2-4], realizing related ideas [5-8] in a concise framework. Space-time is described 
in terms of a 3-1-1- dimensional noncommutative (NC) brane solution embedded in M^°, 
i.e. a quantized space Aie C M}^ which carries a non-degenerate Poisson tensor 9^^{x). 
All matter and gauge fields live on this space-time brane, and there are no physical fields 
propagating in the ambient 10-dimensional spac^ An effective dynamical metric G^'^{x) 
arises on this space-time brane, which governs the kinetic term of all fields more-or-less 
as in GR. This metric is composed of the embedding metric and the Poisson tensor, 
G'^'^ ~ 6'^'^ 9'^^ g^'u'i and defines the effective gravity seen by matter and fields as in GR. 
However, the dynamics of the effective metric is not described by the Einstein equations. 

An essential difference to general relativity is that the metric is not a fundamental 
degree of freedom, but arises effectively in terms of scalar fields describing the embedding 
of space-time M.0 C ]R^°, and the Poisson tensor 6*^^ describing noncommutativity. This 
makes the dynamics of emergent NC gravity somewhat difficult to disentangle. In this 
paper, we obtain approximate solutions which correspond to static and somewhat localized 
matter distributions, having in mind galaxies and their stars inside. The bottom line is 
that large matter clusters such as galaxies are embedded in "gravity bags" or halos, which 
are deformations of the embedding M.e C. M}^ with very long (galactic or cosmological) 
wavelength. These gravity bags turn out to screen gravity at large distances, and enclose 
an effective positive "vacuum energy". Localized matter distributions such as stars then 
induce local deformation of this large-scale embedding, which leads to Newtonian gravity 
within these gravity bags at shorter scales. 

^unlike in other branewoiid scenarios such as [9, 10] 
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The effective vacuum energy inside the gravity bags can be considerably larger than 
the currently preferred value in the ACDM model, which - along with other contributions 
- leads to a significant enhancement of the (galactic) rotation velocities at large distances. 
This might provide an explanation of the observed galactic rotation curves without invoking 
large amounts of dark matter. Moreover, the effective Newton constant is determined by 
the large-scale deformation resp. the gravitational background, and therefore can vary 
somewhat in different locations of the universe. 

The screening of gravity and the different physics of vacuum energy naturally leads to a 
consistent cosmological picture, where the localized matter distributions are embedded in 
a Milne- like cosmology [1 1] . This is in remarkably good agreement with basic observations, 
requiring considerably less fine-tuning than in the standard model. We briefly recall this 



cosmological solution in section 4.4, and show how the present results naturally fit into this 
cosmological context. 

Due to the nonlinear nature of the problem, the results of this paper are somewhat 
incomplete and preliminary. However, the basic results concerning the Newtonian limit and 
the long-distance deviations from Newtonian gravity and GR are expected to be reliable. 
At short distances the approximations are less reliable, and a more complete analysis is 
required before seriously addressing e.g. the solar system precision tests. But in any case, 
it is remarkable how naturally emergent NC gravity seems to reproduce most of the basic 
observations, with far less fine-tuning than in the ACDM model. Given the rigidity of the 
model and its perspective to define a full quantum theory of fundamental interactions, it 
certainly deserves to be studied very thoroughly. 



2 The Matrix Model 

We consider the following type of Yang-Mills matrix models 

Sym = -AtT^T (^^[X'',X'][X''' ,X%aa'Vw + ^W^4X^M/]^ (1) 

where rjaa' = diag(— 1, 1, 1); the Euclidean version of the model is obtained by re- 
placing rjaa' with 6aa'- The degrees of freedom of this model are hermitianj^ matrices 
X^ G Mat{oo, C) for a = 0,1,2, D — 1, as well as Grassmann-valued matrices \Ef which 
are spinors of SO{D — 1, 1) resp. SO{D). The Ta generate the Clifford algebra in D dimen- 
sions. We introduced also an energy scale Aq which gives the matrices X" the dimension 
of length. The action is invariant under the fundamental gauge symmetry 

X'^ ^ f/"^X^f/, ^^U^^W U eU{H) (2) 

^in the Minkowski case, we will assume that the time-hke matrices are anti-hermitian, see below. This 
is consistent with a real metric and action, and will be addressed in more detail elsewhere. 
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as well as a global ISO{D — 1, 1) resp. ISO{D) symmetry, where translations act as 
X"" X"- + c"l. However there is no space-time or geometry to start with; space and 
geometry "emerge" only as solutions or backgrounds of the model. The models can be 
obtained as dimensional reduction of large- super- Yang-Mills theory to dimension zero. 
The IKKT model with D = 10 is singled out by an extended matrix supersymmetry [1]. 

It is easy to see how space (time) arises in such a modej^ Dropping the fermionic terms 
for now, the equations of motion 

[X^[X^X»']]r7,,, = (3) 

admit in particular 4-dimensional noncommutative or quantum spaces M.0 C M}^ as solu- 
tion. This means that we can split the set of matrices as 

X'^ = (X^0O, /i = 0,...,3, z = 4,...,9 (4) 

where the 4 generators X^ are assumed to generate the full matrix algebra Mat(oo,C), 
which is interpreted as space of (noncommutative) functions on M., i.e. Mat(oo,C) = 
Ce{M.)- This is the basic idea of noncommutative geometry. The "scalar fields" 0* = 0*(X^) 
are then functions of X'^. The prototype of such a solution is the Moyal-Weyl quantum 
plane where = iO^'^l, = 0, but we will focus on the case of nontrivial 0* here. 

The basic hypothesis in is that space-time is described by such a quantum space solution 
of We focus on the semi-classical limit of such quantum spaces, indicated by ~. Then 
the 0*(a;) define the embedding of a 4-dimensional submanifolc^ c M^°, and 

[X^X'^]~zr^(x), /i,z. = l,...,4 (5) 

can be interpreted as Poisson structure on A^. In particular, the matrices ~ are 
interpreted as quantization of coordinate functions on Ai. Thus the matrix model provides 
preferred coordinates x'^, which however have no physical meaning whatsoever. From the 
point of view of GR, they essentially "fix the gauge" , disposing of diffeomorphism invariance 
which does not make sense in the matrix model. Since gauge-dependent objects are always 
unphysical, this has no implications on the physical content of the model. 

All physical fields in this model arise from fluctuations in the matrix model around 
such a background (leading to nonabeliaiij^ gauge fields and scalars) and from the fermionic 
matrices Since Mat(oo,C) = Cg{Ai) by assumption, it follows that they all live only 
on the brane Ai, and there is no physical higher- dimensional "bulk" which could carry 
any propagating degrees of freedom. This does not exclude the existence of compactified 
physical extra dimensions in the matrix model, but these are different backgrounds which 
we will leave aside for simplicity here. 

■^This basic observation has been made by many authors including [12-15], and it is obvious from the 
point of view of NC gauge theory [16-18]. 

^cf. [19] for a different approach to NC submanifolds. 

^In the nonabehan case the background solution is generahzed to include a su(n) factor, see e.g. [4,20] 
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Emergent geometry. The Poisson tensor 9^^^{x) not only governs the noncommutative 
structure of TVI, it also plays a crucial but implicit role in the low-energy effective action 
and the metric on M.. We need to assume that it is non-degenerate, so that its inverse 

0~l{x) (6) 

defines a symplectic 2-form on M.. Then the trace on Ma)f:(oo,C) is given semi-classically 
by the volume of this symplectic form, 



(27r)2 Tr/ ~ j d^x p{x) f 

P{x) = {dete-jy/'. (7) 

We can now extract the semi-classical limit of the matrix model and its physical meaning. 
To understand the effective geometry of Ai'^, consider a test-particle on A^^, modeled by 
a scalar field ip for simplicity (this could be e.g. an su{k) component of 0*). In order to 
preserve gauge invariance, the kinetic term must have the form 

S[^] ^ -^r[X^y;][X^¥;]r/,, = -Tr([X^y;][X^v;]r7^.+ [0^¥;][0^¥;]5,,)• (8) 
Expressing the 0* in terms of and using 

[<P\ fiX^] ~ z0^''d,<P%f (9) 

this kinetic term can be cast into covariant form 

S[f] ~ ^ j d^x\G^,\"^G^"{x)d,ipd,ip, (10) 

where [3] 

G^-(x) = e-'' e^^'\x)e'"'' {x)g^,,,{x) (11) 

9f,uix) = df.x'^d^x'' = 7]^^{x) + df,(l)'df,(l)^6,j. (12) 

Here gp.u{x) is the metric induced on A^^ C M}^ via pull-back of rjab- Therefore the kinetic 
term for on A^^ is governed by the metric G^!^{x), which depends on the Poisson tensor 6^^" 
and the embedding metric gfj,u{x). It turns out that the same metric also governs nonabelian 
gauge fields [4] and fermions [3] in the matrix model, hence G^j, must be interpreted as 
gravitational metric. There is no need and no room for invoking any "principles" . 

We note that 

\G^u{x)\ = \g^u{x)\, (13) 

which means that the Poisson tensor 6"^*^ does not enter the Riemannian volume at all. 
This is important for stabilizing flat space, as we will see. Note also that the matrix model 
action ([T]) can be written in the semi-classical limit as 

Sym = -AtTr^-[X%X'][X^',X%aa'Vw ~ / d'xAtp{xUx), (14) 
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where 

Vix) = ^e''G^%x)g,,{x). (15) 

Covariant equations of motion. As shown in [3], the basic matrix equations of motion 
(|3| can now be cast into a covariant form as follows: 

= (16) 
□gx'^ = (17) 
V^(e'^^-;) = e-^'G^^e^^d.vix) (18) 

where we consider x'^ ~ as a scalar function on A4 , consistent with the ambiguity of the 
splitting X"" = {X'^, 0*) into coordinates and scalar fields. Here Vg denotes the Levi-Civita 
connection with respect to the effective metric C^". In particular, such on-shell geometries 
imply [3] 

= (19) 
for the preferred matrix coordinates x'^, which in general relativity would be interpreted 



as gauge condition. Furthermore, it turns out that (18) is in fact a consequence of a ma- 
trix Noether theorem due to the translational symmetry X" —>■ X"- + c"ll, and is therefore 
protected from quantum corrections [4] . It provides the relation between the noncommuta- 



tivity 9^'^{x) and the metric G'^'^. Since (18) has essentially the form of covariant Maxwell 
equations coupled to an external current, it should have a unique solution for a given 
"boundary condition" 

0^v{x) = const for |x| ^ oo (20) 

up to radiational contributions, which play the role of gravitational waves here [2,5]. 



Relation with string theory. The IKKT matrix model was proposed originally as a 
non-perturbative definition of IIB string theory on M}^. From this point of view, Aig could 
be interpreted as a brane with open string metric G'^", while Qf^i, could be viewed as closed 
string metric in the lOD bulk. Indeed there are also other solutions of the matrix model, in 
particular lOD solutions. Notably graviton scattering has been studied in this and related 
matrix models; for an incomplete list of references see e.g. [1,12-15,22-25] and references 
therein. Most of this work was based on a different kind of block-matrix backgrounds, 
and to simple geometries in the NC case. The essential point here is to consider general 
4-dimensional NC brane solutions, without physical lOD bulk. On such a background, the 
matrix model can be viewed as NC gauge theory [16,17], which includes gravity as we 
have seen. In this way the strength of string theory (notably the good behavior under 
quantization) is preserved while the main problems (lack of predictivity) are avoided. In 
particular, the matrix model should be viewed as background-independent here, since there 
is no physical space-time to start with. 



6 



A natural question arises how the present resuhs relate to the BFSS matrix model [23] 
of (M)atrix theory, which involves 9 matrices depending on a classical "time" parameter. 
This model is not unrelated to the IKKT model as pointed out in [1], and some version of 
the present mechanism should be realized there as well; this should be studied elsewhere. 
We only want to stress here that the space-time brane Ai must be entirely noncommutative 
in the present framework, i.e. 6^'^ must be non-degenerate. Nevertheless a conventional 
physical picture emerges at scales larger than Ajvc and there is no need for a classical 
time parameter. Hence the IKKT model is more natural for the mechanism of emergent 
gravity. It is also very interesting that some evidence for 4-dimensional space-time (rather 
than some other dimension) emerging in that model has been found [26]. 

The topic of membranes and matrix models has of course a long history, cf. also [27-29]. 
It is also interesting to compare the present approach with other models of emergent gravity, 
see e.g. [30] and references therein. 



2.1 Self-dual solutions and (/^^j^ = G^j, 



A simple class of solutions of ( 18 ) is given by "self-dual" 2-forms 6^^ (x) which satisfy 



0, 



(21) 



To see this, consider (at a point x) a local coordinate system where g^i, = diag(s, 1, 1, 1) 
(with s = ±1 in the Euclidean resp. Minkowski case) and 6'^'^ has the form 












— y/sa 




































(22) 



this can always be achieved using a SO{A) resp. 50(1,3) transformatioi]|^ at the point x. 
Clearly this is (anti-) self-dual -kO'^ = ±^/s9~^ if and only if = 1, where * denotes the 
Hodge star and 6^^ = 9~ldx^ A dx'^ . On the other hand. 



2-2 2 



diag(sa^, a 
. Therefore -kO^^ 



at the point x, hence e~'^7] 

dfMu = Gfj,^, or equivalently e'^ = r]. It is then easy to see that in this self-dual case ( |18[ ) 
reduces to V^6~l = 0, which holds identically since d-k = d9~ 



(23) 

±^/s6~^ if and only if 
is 
0. 



X 



X 



Such (anti-)self-dual closed 2-forms 6 ^ with constant asymptotics 6 
oo always exist, at least on asymptotically flat spaces. This can be understood by 



6^^ as 



^Note that we assume that is imaginary in the Minkowski case. This might appear strange, but 
it is natural having in mind a Wick rotation x*^ = it, and it is essential for g^j^ = G^i^. This should be 
addressed in more detail elsewhere. 
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interpreting as sourceless electromagnetic field with constant field strength at infinity. 
Indeed, we only have to solve d-k F = 0, F = dA with suitable asymptotics F ^ F as 
r oo, and define 6'^ to be the ( ant i-) self dual component of F. 



In this self-dual case g^u = G^^, the bare matrix model action (14) becomes 



Smm = Atj d'x pr] = A.^ j d^x pe'' = j d^x ^J\G^\K = j d'x ^/\^\At (24) 

which is precisely the form of the vacuum energy resp. brane tension, interpreted as 
cosmological constant in GR. We collect all such (bare and induced) terms in 



S,ac = -2 J d''x^\g^,\K\ (25) 
denoting with the sum of Aq and the induced quantum-mechanical vacuum energy; the 



above sign is essential for stability reasons, cf. section 4^ We also recall that the effective 
Yang-Mills coupling constant for the basic SU (n) gauge fields on A^^ is given by [4] 

^^^, = Ao^e--~Ao-U^^^, (26) 

so that qym = 0(1) amounts to 

Ao ~ Ajvc- (27) 
Let us briefly discuss deviations from the self-dual case g^^ ^ G^^. These arise from 
(abelian) variations Sx'^ = 9^^ Ai, of the tangential matrix degrees of freedom resp. "would- 
be U{1) gauge fields", which lead to 56~l = F^iy where F^i. is the corresponding field 
strength. The corresponding propagating degrees of freedom behave as gravitational waves, 
with metric fluctuations h^^, = Gy^iO^ ''Fp^ + G^^S^ ^Fp^ + \GpyFpr^6'^ [2,5]. In the presence 



of matter, we expect additional source terms on the rhs of equations (16) - (18). While 
matter is not charged under this gravitational U{1), it is expected to induce e.g. dipole 
(and higher multipole) excitations of Fp^. This would lead to short-distance modifications 
of the metric which decay more rapidly than -. We therefore simply ignore such effects 
here, and consider only the simplest "self-dual" solutions with gp^ = Gpy. 



2.2 Static metric deformations 

We focus on deformations of the embedding C MP through the scalar fields 0* with 



embedding metric ( 12 ) 



gi,u = Vf^u + dp(l)'du(p' 

= Vi^u + hpy, (28) 

and assume that Gpy = gpu, i.e. that 9'^'^ is self-dual w.r.t. gpu as discussed above. This 
is an ansat:^ which may not be completely appropriate, but it is expected to give the 

^In fact one of the main difficulties in this context is to understand and disentangle the effects of the 
embedding versus the NC structure 9'^'^ . We hope that this paper helps to clarify this issue. In particular, 
the ansatz in [2] with trivial embedding but nontrivial 6'^'^ may not be realized in the presence of matter. 
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correct long-distance physics in the static case since corrections due to G^i, ^ g^^y lead to 
be short-distance effects and gravitational waves. 

We furthermore focus on static metrics g^^, corresponding to static and somewhat 
localized matter distributions. One goal is to obtain the analog of the Schwarzschild solution 
at least in some regime. Since hfj_u is quadratic in 0, we cannot restrict ourselves to linearized 
fluctuations in around flat Minkowski space. We must include at least quadratic terms 
in (p, which makes the analysis somewhat non-trivial. Nevertheless we will find a class 
of solutions which appear to give the appropriate description of objects such as galaxies 
sparsely embedded in flat Minkowski space. This also fits naturally in the context of the 



cosmological solutions [11], as discussed in section 4.4 



Some insight can be gained by observing that the metric fluctuation h^j^^ is closely related 
to the energy-momentum tensor of the massless scalar flelds 0*, 

^/t = df,(f)dy^ - ^Vf^uiv^'^da^dp^) = V - ^r]^y h, 
h = t^Kp = -r^^''Tf:^. (29) 

In particular, the effective Newtonian potential is related to the energy density of a free 
scalar fleld, 

2U{x) = -/loo = -Tot, + \h^ -^Tot, (30) 

assuming Tf- in the last step. Hence a static metric fluctuation corresponds to a a 
scalar fleld excitation with static energy-momentum tensor. Furthermore, since in vacuum 
satisfles the free wave equation ( [T6| , it will decay like ^ outside of matter distributions, 
suggesting a U{r) ~ ^ behavior for the gravitational potential. This appears very bad at 
flrst, since this is not the usual f/(r) ~ ^ law of Newtonian gravity. 

The resolution of this puzzle is as follows: the relevant non-smg"u/ar harmonic excitations 
necessarily have some flnite wavelength, 0* ~ EBl^gj^^t^ These are localized excitations 
resp. "gravity bags", which we will argue to have very long (astronomical) wavelength 
uj. This leads to a long-distance screening of gravity with U{r) ~ ^. However the crucial 
point is that inside these gravity bags (in particular inside of galaxies, say), the local matter 
distribution does indeed lead to the standard f/(r) ~ ^, as we will show. Newtonian gravity 
therefore arises as a "short-distance" effect on the harmonically embedded spacetime brane. 



Basic "rotating" embedding. We consider brane embedding C which are 
small deformations of flat Minkowski space, with the following structure: 



X 



A 



/ x'] 








) 





( 



t 

COs(co't) 

sin (cut) 



(31) 
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where g{x) is independent of t and g{x) ^ as r ^ oo. There might be additional 
components in this ansatz. We stress that the 0* comprise two space-like components 
which rotate with t, rather than just of its component^ We will see in section 4.4 that 



due to its flat asymptotics, this type of embedding fits naturally into the context of the 
cosmological solutions obtained in [11]. An important feature of this ansatz is 

docp'dicP" = 0, (32) 

leading to a static metric 

cis^ = -{l-uj^ g^) dt^ + {5ij + digdjg)dx'dx^ (33) 

This is consistent with the fact that there is no energy flux associated to standing waves. 
There might be several components of g —>■ g^ whose contributions will simply add up; this 
will not alter the essential conclusions below. Clearly the rotating embedding with uj 
must be responsible for Newtonian gravity, which cannot be reproduced with a purely static 
embedding. Additional components might be necessary to comply e.g. with the detailed 
solar system constraints. The embedding functions 0* must be determined by solving the 
equations of motion in the presence of matter, which will be done below. Note that the 
preferred "matrix coordinates" x^ automatically satisfy the e.o.m. OgX'^ = in vacuum 



(16), (17); therefore we will only consider the equations of motion for 0' below. 



3 Coupling to matter and effective gravity 



3.1 Effective action and equations of motion with matter 



In this section we derive the equations of motion for gravity coupled to matter in the 
semi-classical limit, assuming g^^^, = G^^, as discussed above. Our starting point is the 
semi-classical effective action of the matrix model (25) together with the action for matter, 
and the Einstein-Hilbert action R[G] = R[g] which is induced at one-loop[^ 



S = j d'^x^/lil {AlR - 2Af) + S, 



matter- 



Now recall 



9\ 

5 I x/\g\R 

'^'S'matter 



6^-detg = ^^\g''''6g^ 
-VWlG^'Sg,, 



flUI 



(34) 



(35) 



^Thc time-like component (s) could have similar nontrivial embeddings. This should be elaborated 
elsewhere. 

9 Assuming a cutoff A4, which is the cutoff for = 4 SUSY in the IKKT model, see [31, 32] 
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where Q'^" = R^^ — \g'^^R is the Einstein tensor, and T^'^ is the energy-momentum tensor 
for matter (recall that matter and fields couple to the effective metric essentially in the 
standard way). The crucial point is now that the fundamental geometrical degrees of 
freedom are not the g^i^, but the embedding fields 0* as well as resp. 9^'^ . In the self- 
dual case g^j_v = G^^, the most general variation can thus be decomposed into variations of 
dfiu = Vfif + d^cjfdyctf = G^u and variations of 9^'^. The e.o.m. for 9'^'^ are satisfiec^" due 



to (21). The variation of S with respect to the fundamental fields 0* can be written using 



Sg,^u = df,(f)'dj(f)'' + 9^50*9^0* as 

6S = j d^x^\5g^,W = -2j 6<p%i^/\i\'^^n^u(p' (36) 
up to boundary terms, where 

n^"^ = SnT^"" - Alg^" - A^g^"". (37) 
This leads to the equations of motion for 0* 

^^,{^M'>^^''^A) = 0, (38) 
which using the identity Vf.V = -^d^{^V^) can be written as 

A{Ugct> = [SnT^" -Alg^'')V^d,<P + 87i{V^T^ndu<P (39) 
recalling V^Q'^'' = 0. This equation has 2 types ("branches") of solutions: 

1. "Einstein branch" : 

Clearly every solution of the Einstein equations 

Ag^"" + 0^" = 87tGT^"' (40) 



is also a solution of H^'^ = (39), upon identifying 



^ = § = At (41) 

Indeed using general embedding theorems [33], any solution of the E-H equations can 
be realized (locally) using embeddings in D > 10. This should provide an interesting 
realization of Einstein gravity within matrix models, with technical advantages over 
other approaches in particular for quantization. Recall however that we assumed 
Gfj^i, = g^y, which will not always hold; thus there will be some modifications of the 
Einstein equations due to 9^'^{x). On the other hand, this Einstein branch would 
presumably lead to the same problematic aspects of GR, notably fine-tuning issues 
in cosmology. In this paper we will focus on the second kind of solution, which leads 
to significant and very interesting deviations from GR at long distances. 

^"This however amounts to neglecting matter contributions to as discussed before. For geometries 
with non-selfdual 9^^^ this derivation must be refined. 
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2. "Harmonic branch" : 



Besides to the above solutions, there are additional solutions of (39) with 

d^i^M^^'d^'P) = 0, ^^"^ + (42) 
The prototype of such a solution without matter is flat Minkowski space 0* = with 



Ai > 0. More generally if the vacuum energy dominates the matter density, then (39) 
reduces to ~ 0, leading essentially to minimal surfaces which will be deformed in 
the presence of matter. This in turn leads to the near-realistic cosmological solutions 
of FRW type found in [11] with Milne-like late-time behavior, which are stable and 
largely insensitive to the detailed matter content. Another attractive feature of this 
harmonic branch is that its quantization should be comparably simple, as the em- 
bedding fields 0* are governed in vacuum by a simple action with positive excitation 



spectrum, cf. section 4.3 



We will simply assume in the following that lS\Q'^^ can be neglected in the static case. 
This is however far from trivial, since the simple ansatz (31 ) will generally not lead to Q'^^ ~ 



0. There are 2 ways how this might nevertheless be justified: either a more sophisticated 
ansatz taking into account Q^^ will imply Q"^^ ^ 0, or A^ is much smaller than the Planck 
scale. The latter is in fact very appealing as we will see. In any case we focus on the 
harmonic branch in this paper, since it leads to very interesting long-distance modifications 
of Newtonian gravity which is the main topic of this paper. 



3.2 Perturbations of flat Minkowski space-time 

For simplicity we focus on perturbations of flat Minkowski space, and only consider the 
case of "weak" gravity, i.e. 

9i,v = VfMiy + h^,y (43) 



keeping only terms linear in h^i, = 9^0*9,^0*. The e.o.m. (39) then simplifies to 



usmg 
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= ^ T^"9^9,0, T^^ = T^" - ^g^" (44) 

V^T^" = 0. We also replace Dg ^ to leading approximation. In the case of a 
static mass distribution p, we assume furthermore that 

Too = Pix) = p"^'^"^''(x) - p^"™(x) > 0, p,ur. = ^g"' (45) 

^^The conservation law for T^'' is not entirely evident in the matrix model, in particular for fermions 
which have a non-standard spin connection [34] . Nevertheless it is expected to hold at least to a very good 
approximation since the matter action almost coincides with the usual one. 
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and Tij as discussed abov^^ Then (44) can be written using □ 



-d? + A as 



in 



A0 = (1 + ^pix))d'o<P- 



Note also that 



(46) 



(47) 



so that /i^j, satisfies the "harmonic gauge" d^h^y = \dyh up to corrections due to matter 
and curvature. 



3.2.1 Harmonic gravity bags 

We are mainly interested in static spacetime-geometries in this paper. Thus consider the 



following localized excitation of the embedding (31) 



nx,t) = g{x)e^-' = g{x)h2^^^^^^ ) (4J 



with some very small tu, leading to the static metric (33) 

ds^ = -(1 - ^ 1^^,, ^ digdjg)dx'dx^. (49) 



For p = and neglecting the curvature corrections AlQ^^, the equation of motion (46) 
becomes 

□0o(a;) = 0, Agoix) = -uj^go{x). (50) 
Consider first the case of a spherical wave, where this equation reduces to 

dr{r^g') +ujh-^g{r) =0. (51) 

The unique spherically symmetric solution which is regular at the origin and decays as 
r ^ oo is given by 

sin(ci;r) 

9o{r) = go , 

ur 

M^) = S.<^)^'=<)'>^-)i'^^) (52) 



with radial wavelength given by 

J^uj — ~- 
00 



(53) 



^^One might think that due to the Einstein equations, p = Pmatter ~ Pgrav = would drop out of the 
harmonic equation. This is not the case here, as the Einstein equations do not hold in this form. 
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Near the origin, we can write 

9oir) = (70^^ ~ ^70 1 - + ... , 

ior o 

g'{r) = go( ^ ^ go{~-uj\ + —u'r^ + 0{r')). 54 

Note that g' is regular at the origin, and decays hke - as r ^ oo. The effective metric (33) 
is 

ds^ = -{\-u? g{rf) de + (1 + {g'f)dr^ + r^dn\ (55) 

which allows to read off the effective gravitational potential: a static test particle in this 
metric perturbation feels an effective gravitational potential 

^oo = -(l + 2f/o), 2[/o(r) = -u'g{rf = -uj'gl('^^)\ (56) 



which satisfies 



Uo{r) ~ r ^ oo (57) 



r 



AUoir) = gluo\ r~0 (58) 

The last equation will be interpreted later in terms of an effective gravitational con- 
stant/vacuum energy Aeg (72) inside the bag. For large distances, U{r) leads to a rapidly 
decreasing attractive gravitational force with range L^^. This might appear strange at first, 
since such vacuum excitations do not exist in general relativity. However, they are not 
expected to survive in this naked form: Due to their attractive gravitational force, mat- 
ter will tend to accumulate inside, and these "gravity bags" will typically enclose matter 
with p 7^ 0. In particular, large clusters of matter such as galaxies will be embedded in 
such "gravity bags" . The essential point is that the matter within such a gravity bag will 
experience Newtonian gravity, as we will show next. 



3.2.2 Spherically symmetric mass distribution and Newtonian limit 

We now show how Newtonian gravity arises in the simplest case of spherically symmetric 
localized mass distributions, due to a local deformation of the above harmonic gravity bag. 
This will lead moreover to significant long-distance deviation from Newtonian gravity and 
GR. Basically, gravity is confined within the bag, which also contains an effective vacuum 
energy. The case of a more general mass distribution will be studied in section |4j This will 
clarify the significance of the gravity bag which we simply assume here. 

Consider a spherically symmetric mass distribution around the origin within the radius 



tm- For r > rjvf, the corresponding solution of (46) with the "static" ansatz (48) must have 
the form 



cj)' = g{r)e 

, , sin (cur + 5) , , ^, sin(5), 

9{r) = 90— ~ go{cos{6) + 59 
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assuming cot <^ 1 in the last expression. This phase shift 5 7^ is the key for obtaining 
Newtonian gravity. It leads to the effective metric 



5-00 = 



l2~ 



= -(1 - gy cos(25)) + ^ ^ + --^^^^^ - gy sin(25) 

T V O 

-\gycos{26y + 0{r^), (60) 

1 I '2 
grr = l + g 

1 5fnu; sin(25) g'^sin^iS) 5fnSin^(5) 2 o ^ • /^e-x 

= 1 + -^ — + ° 2 + 24 + T^^oVsm 2(5 r 

+ ^^^^^cos(2(5)r2 + 0(r^). (61) 
9 

This corresponds to a gravitational potential f/(r) = —^uj'^g^, which for intermediate dis- 
tances 

sin5 < cjr < 1 (62) 
is well approximated by a ^ potential with a constant shift, 

TTf \ 1 2 2_1 2 2 /ox^ ^o^sin(25) 

U{r) = --L0 g ^ -g^uj cos{2S) — . (63) 

Thus all we need to obtain Newtonian gravity is 5 ~ M, which is very intuitive and indeed 
correct as shown below. 



Boundary condition. Now consider the region near the origin where p(r) > 0. Equation 
(46) gives 

= c.V(l + ^p)<7(r) + 9,(rV). (64) 
We focus on the region near the origin, wher^^ g ^ g{0). Then 

r^g'ir) = -u'^ I dr'g{r)r'^{l + ^p) 

-u'giO) I dr'r'\l + ^p) 
2M(r)~ 



-^^^(0) ( V + ^ ) (65) 



hence 

.'W /. 2Mg) . ^^^^ 



^(0) V3 r^Af 

^■^this will be justified in section |4j noting that is due to a large-scale background structure which 
dominates the local perturbation due to M. 
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where M(r) is the mass inside the sphere with radius r (including possibly curvature 
contributions). Assuming that p{r) = for r > tm, we can match this with (59) which is 
valid for r > tm'- 

, . sin(co'r + 6) 

9{r) = 90 

r 

q' , ^, , cos(a;r) sin(co'r), ,cos(ti;r) sin(cijr), 

— = cos 5 ^ — ^ -sm 5 V^ + cj — ^ — ^ 

go J. ^ 

^ ' [1 + —-) - -u^ cos{6) r + 0{uj^r^) (67) 



r2 ' 2 ' 3 

assuming cur <^ 1. Note that the singular terms are misleading, since this expression is 
valid only for r > vm- Combining this with (66) gives the matching condition 

g'{r) sin((5) 1 ^ , 2^(0) fr 2M 

— ^ ^^ + -^'cos 5 r = uj^^l + 

where we neglect the constant term since ur <^ 1. This implies the two conditions 

giO)2M 3 . 

-uj = sm{d), 



90 

9{0) 



-- cos((5). (6^ 

90 



Thus the time-component of the effective metric (60) is 

1 + 2t/o - 2^ + ^MG.;V - ^A.^r^ + {—f^ + 0(r3) (69) 
r 3 3 r 2t/o / 

where we define 

Uo = -\9lu^\ (70) 

G - ^--4f^o^ (71) 

Aeff = -]^GK\ = 2Uouj\ (72) 

(note that G is naturally small given that is large), assuming 5 <C 1 and cur <C 1. 
Let us discuss the importance of the various terms. The ©(M^/r^) term 

9l^^n^_MG 1 

r^ r ^ 2f/n' ^^^^ 
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can be neglected compared with the Newtonian term provided 



sin (5 < ur i.e. < glu"^ = \2Uo\, (74) 

which we assume for simphcity. This means that the Newtonian potential due to the 
local mass M should be smaller than the background potential Uq due to the harmonic 
baj^ This also implies that the vacuum energy term (as well as the Newtonian potential) 
dominates the linear term, 

lAeffr^l = 2|f/o|cuV » MGcuV = — (75) 



which we will omit henceforth. Thus the time-component qqq of the effective metric (69) 
has approximately the form of a Schwarzschild-de Sitter metric with a constant shift [35], 



900 



;i + 2t/o - ^ - \Ksr') (76) 



assuming < r < and dropping the linear term l^^cu^r^. The Newtonian term 
dominates the Aefrr^ term provided the vacuum energy 

EUr) ■■= « M, (77) 



(cf. (119)) inside r is smaller than the mass M. We then obtain Newtonian gravity with 
potential 

Uir) ^ f/o-— . (78) 
r 

This will be generalized to the case of an arbitrary mass distribution in section |4j Finally, 
we note that 

^° - "16;^ ^^^^ 

can be interpreted as Newtonian potential due to the vacuum energy contained within the 
harmonic bag of size L^. 

The basic result is that Newtonian gravity arises at intermediate scales, with important 
long-distance modifications. Notice that the precise form of the induced gravitational action 
was never used up to now, rather gravity arises through a deformation of the harmonic 
embedding which couples to T^'^ . Thus the mechanism is quite different from GR. 



^'^Having in mind e.g. a star within a galaxy; this will become more obvious in section [4] in the context 
of a general mass distribution. 
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3.2.3 Deviations from Newtonian gravity 



We have shown that under the above assumptions, goo has essentially the form of a 



Schwarzschild-de Sitter metric with an apparent negative cosmological constant Aeff (72) 



which is related to the vacuum energy. However, note that the sign of AgfF is different from 



GR (41 ). This is not a mistake, but underscores the fact that the physics of vacuum energy 
is different here from GR. Inside the gravity bag, the vacuum energy > contributes a 



positive energy density to the gravitational potential. Note also that (|58|) for the harmonic 
gravity bag can now be written as 
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Af/o(r) = -Aeff = 4716-^ (80) 

87r 

corresponding to an effective energy density due to A^ However for very large distances 
> -^cj; the effective metric approaches the harmonic behavior 

, In ,sin^(co'r) , , 

U{r) ~ --uj^gl^^ (81) 

which is rapidly decaying and oscillating, smoothly merging with the flat large-scale metric 



9fj.u{x) — > resp. the Milne-like cosmology [11] as discussed in section 4.4, In particular, 
we will see that cosmology does not lead to the usual stringent constraints on the vacuum 
energy. However we obtain an upper bound for Ai e.g. due to solar system constraints (for 
example Ai = 0{eV) would work, but certainly not Ai = 0{TeV)). 



A plot of the full U{x) = —i^oj g^x) in comparison with the terms in (76) is given in 



figure [T]for uj = 0.1, go = 1,6 = 0.1 hence Uq = —0.005. This clearly shows the dominating 
Newtonian form for small r. The vacuum energy term takes over for larger r, until the 
potential is cut off effectively at L^. Note that the amplitude of these oscillations is very 
small and thus may be very hard to detect, although in principle such oscillations should 
be present at very large distance from isolated spherical galaxies. They might be masked 
for various reasons such as matter distributions or the superposition of slightly different Ui. 



Now consider the radial part of the effective metric (61), which can be written as 



12GM 1. . CM 1 GM , . 



There is a strange factor | in (82) compared with general relativity. However, one should 
keep in mind that this metric is appropriate only for isolated masses resp. galaxies, but not 
for small perturbations within galaxies such as the solar system. The latter will be studied 
in section |4| confirming the factor | for stars within galaxies, while the more singular terms 
are smaller in that case. This may be a challenge for the solar system constraints. However, 



^^recall that in GR, a negative A leads to an additional attractive gravitational field as above, due to 
AfJ = 47rGpniattcr " A, See [35]. 
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U(r) 




Figure 1: Gravitational potential U{x) compared with the Schwarzschild-de Sitter potential 
Uq — — |Ar^ (dashed line) and long-distance oscillations, for u = 0.1, = 1,S = 0.1. 



there will be short-distance corrections e.g. due to 9^'^{x), and a more complete analysis is 
required before a reliable judgment can be given. 

Note incidentally that the metric is regular but becomes degenerate as qqq — » 0. How- 



ever, then the approximation of linearized gravity in (44) is no longer valid, and a more 
careful treatment is required. 



(Galactic) rotation curves. Now consider (non-relativistic) orbital velocities. For small 
distances r < L^, it is given by 



v = VlN= W2 — (1 + -jjj^) - ^Kfjr'. (83) 

This decreases like r^^^"^ as in Newtonian gravity as long as E^acir) < M, but for Eyaci^) ~ 
AttM it starts to increase linearly like v ~ A/|Ae//| r until r ^ L^; recall that A^jj < 0. At 
that scale, the velocities decrease again with oscillating behavior and 

/ \/2 

V = VlPr ^ —^0 (84) 
r 

A plot of this orbital velocity for the exact U{x) and the Newtonian approximation is 
given in figure |2| for a; = 0.1, = 1;^ = 0.1 thus Uq = —0.005. Note that these simple 
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Figure 2: Orbital velocity v{x) for a central point mass compared with Newtonian case: 
Newtonian domain, enhancement and cutoff, for uj = 0.1, go = 1,6 = 0.1. 

relations hold only outside of the mass distribution, and will be modified by the presence 
of mass e.g. in the halos of galaxies, and by the onset of the harmonic long-distance decay 



(84). The result is similar to that of an Schwarzschild-de Sitter geometry with Aeg given 



by (72), combined with a harmonic screening at r ^ L^^. It remains to be seen if this 
indeed allows to explain the observed galactic rotation curves; however qualitatively, the 
above behavior certainly goes in the right direction. What is particularly striking is that 
it naturally predicts a slightly increasing rotation curve, which is indeed often observed. 
This may be masked by other effects, in particular the non-trivial matter distribution. 
In fact, the possibility that the galactic rotation curves might be explained in terms of a 
cosmological constant has been proposed in the literature [36]. It was argued to be feasible, 
provided A^s can depend on the individual galaxy (the value Aeg — 5 x lO^^^cm^^ was 
given as a typical scale). Of course this does not make sense in conventional GR, and it 
is inconsistent with cosmological constraints in the ACDM model. However the latter are 
irrelevant here as discussed below, and G = is dynamical here and may indeed depend 
on the individual galaxy. It would be extremely interesting to study this in more detail. 



Some estimates. We certainly want to preserve 5^00 < 0, which implies that 

1>2.„.,.V.«^^^. ,85) 
This allows to express the vacuum energy as 

A^ = v^C^Apianck = V2U0 (86) 

To get some idea of the scales we consider our galaxy. Assuming ^ -Z^gaiaxy ~ 10^°m (so 
that the oscillations occur roughly at the size of our galaxy), this gives 
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To be specific assume that the potential energy at the center or the galaxy roughly coincides 
with its value given by Newtonian gravity, 



1 QlS™ 

Uo ^ MG/L,^u.y ^ ^ 10-^ (88) 

where M ^ lO^^fcf?. Then we get a bound on Ai of the order of electronvolt. We will see 



in section |4.4| that the cosmology in this model is not very sensitive to the actual value of 
Ai. Moreover, A4 could be much smaller than the Planck scale in this model as discussed 
below. Thus the fine-tuning is considerably milder than in the standard ACDM model 
where Ai ^ 2 x l{)~^eV. 



4 General matter distribution 

In this section, we generalize the above considerations to the case of an arbitrary but still 
somewhat localized mass distribution. The approach is a bit different from the one on the 
previous section, which also provides a confirmation of the above results. 



Consider a large harmonic "gravity bag" solution (52) with size L^^ of the order of a 
galaxy, say, which satisfies 

Agoix) = -u'goix). (89) 
Now add (non-relativistic) matter with p(x) = Tqo ^ and Tij ^ into this bag; indeed 



matter will tend to accumulate in the bag due to its attractive gravitational potential (58). 
If the localized matter density p{x) (consisting of stars, etc.) is not too large, this will lead 
to a deformation 

(j)\x,t) = g{x)e'-\ 

g{x) = go{x) + 6g{x) (90) 

where \6g\ <^ go{x) ~ go{0) is varying on short scales according to p(x), while go{x) is 
slowly varying at the scale L^. Notice that this is precisely the splitting in ([59|), and 



the condition \6g\ ^ go{x) corresponds to (74). Thus go{x) is the embedding due to the 



average, overall mass in the galaxy resp. the largest cosmic structures, and 6g{x) due to 
small local perturbations such as stars, as indicated in figure |3j This dominance of large- 
scale structures is vindicated by the fact that the typical velocities of large-scale structures 



in cosmology are larger than for small-scale structures, cf. (83) 
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Figure 3: sketch of embedding function g{x) with short-scale perturbations and long- 
distance oscillations 



4.1 Metric deformation and Newtonian limit 



In this situation, it makes sense to linearize in Sg{x); this is the crucial step. The corre- 
sponding metric is static and has the form (33) 

ds'^ = -(1 - g^) dt^ + {5ij + digdjg)dx'dx^. (91) 

We want to derive the Poisson equation for the corresponding gravitational potential 
U{x) = —^u^g^. Thus consider 



(A + u;^)U = -u;^g{x){A + u;^)g{x) - d.gd^g. 



(92) 



Now recall the equation of motion (46) 



{A + uj^)g = -^TXUJ^{jg{x) 



Af 



Using (89) this gives for the fluctuations 



(A + u')5g = -S-Ku^-^^igo + Sg) ^ -Snu 



.2_P_ 



go- 



(93) 



(94) 



Now we use the assumption that p is some small mass distribution inside a large background 
bag go which is slowly varying. This means that 



digo 
di6g 



0{rAgo) = 0{ruo'^ 



0{fA5g) = 0{ruj^ 



9o) 



(95) 



where r <^ L^^ denotes the distance from the "center" of the bag, and r is the distance from 
p. Putting this together and using tu^r^ <^ 1, we obtain 



,p(x) 



Af 



- uj^6'W,gd,g 



iiTuj gQ 



?M+o(rW^o^(l + 87r-^)) +0{5g' 



(96) 
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This shows that the terms with first derivatives can be neglected assuming r^iJ^ -C 1, and 
the term 0{5g'^) can be neglected for small p. Identifying the Newton constant as 



G = 2gl 



(97) 



we obtain the Poisson equation 



MJ ^ AuG {p{x) + 



At 



57r 



(9^ 



with an effective vacuum energy which could be interpreted as apparent negative cos- 



mo logical constant Agfj. This confirms the results (71), (72) of the previous section. If the 

(99) 



matter density is much larger than the vacuum energy, 

P 



Af 



> 1 



this gives indeed the usual Poisson equation of Newtonian gravity. 



AU 



A-nGp. 



(100) 



As a check, consider the case of a small localized mass M at the origin. The fluctuations 
^9 9o can be determined exphcitly using (94), which implies (A + u;^)(5(7 ~ —Att^^ 
Outside of the mass distribution this gives 



5g 

9 

GM 

^9l 



9oa 



sin(u;r + 6) 



/sm(Ljr) sin(co'r + 6) 
9o + S9 = 9o{ + a ^ 



\ iur 



90- 



sm{ujr + 6) 



uor 



a sin((5). 



(101) 



This results in a metric which is essentially the same as (60), (61), replacing 6 = a6 for 



small a. However, a more complete treatment including corrections due to 6^'^(x) and 
more general "gravity bags" is required before a detailed comparison with the solar system 
constraints can be performed. 

We conclude that localized matter p inside the gravity bag 0o is subject to (emergent) 
Newtonian gravity, with a dynamically determined gravitational "constant" G given by 



(97). For example, two stars or planets would lead to local perturbations 



9{x) = 9o + S9i{x) + 6g2{x) 



(102) 



where 5gi are perturbations due to object pj. They both see the same go and tu, thus 
the same gravitational constant G, and Newtonian gravity is recovered on scales shorter 
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than L^. However at very long scales and near the border of the galaxy, the effective 
gravitational constant might differ; this would require more detailed modeling. 

Recall that this mechanism is necessarily non-linear and therefore somewhat non-trivial, 
but nevertheless it is quite robust: the essential ingredient is the gravity bag go{x)e^'^* 
which must be slowly rotating and have large amplitude. As discussed before, it seems 
unavoidable that such a bag forms around large matter clusters in the static case. This 
leads to an effective gravity which is at least close to what we observe, and might provide an 
explanation for the rotational curves in typical galaxies without resorting to "dark matter" . 
An obvious question is how the parameters u and qq, and thereby G are determined; we 
will give an argument below why G should be at least approximately the same in different 
galaxies. 



4.2 Ricci tensor and relation with general relativity 

In order to compute the Ricci tensor in the above situation, we assume for simplicity that 
the mass distribution p is in the center of the background gravity bag, so that 

di(f)o ^ diQo = 0, 
didjgo ^ --dijuj'^go (103) 

and go{x) ^ go ~ ^uj'^gor'^ for small distances. The space-like components of the metric 
hij = digdjg can then be written as 

hij = digodjgo + 2digodjSg + 0{5g'^) 
1 2 

= -^gl^^XiXj - -goUJ^Xidj6g + 0{6g^). (104) 



The Ricci tensor can be computed using the on-shell relation ( [47| ) 

d^hxu = \d,h + □0(9,0) = \d,h - Sttcu^ -^0(9,0) (105) 

Zi 1 



hence 



so that 



d^d'-hx, = ]^d^d,h + n<p{d^d,<p) + d^in<p)d,(t) (106) 



R^, = 9,r°, - d^Ti^ = ^ (-□ V - d,d,h + d^d^hx. + d^d^hx,) 

= -(a^a»(a,aA0) + (n0)(a^a,0) 

R = -(9^9V)(5A</') + (n</')' 

= a;V-2a;2aVaA0-a^aWj^7-(87r^'|4)V (107) 
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using (9o0(9o0 = 07^0^ = —0^00 and (103)). This is hard to evaluate in general. We can 



write it as sum of contributions R^i, = B^} + R^,^v + where 
4°^ = -(5m^Vo)(9.9a0o) 

R^^} = -(5^9>o)(a.aA50)-(9,9^0o)(5A'50) + (n^'^)(^A0o) 

Rfu = -{d^d^S<P){d,dxS^) + {DS<P){d^dM) (108) 

since □00 = 0. is the ' 'vacuum" contribution due to the background bag, is the 
desired contribution due to matter linear in p, and R^^} = 0{p^) is expected to be negligible 
for small p. Consider first the "vacuum contribution" R^^^} which applies whenever p = 0: 

Harmonic bag contribution. The contribution of the harmonic bag 0o is given by 

<J = -(5m^>o)(9.9a0o), R = -(9^9>o)(5^5a0o). (109) 

,2x2 _ _Aa2, 



This can be written using 9o09o0 = to (f) = — 09g0 and (103) as 



Roo = -o^29Vo5a0o ^ -^'^7o = -^GAt, 

Rij = -{didk(f)o){djdk(f>o) ^ - l-6iju;^gQ = -^SijCA^ 

y to 

R ^ ^-coY, = ^-GAt. (110) 

Therefore R^^u = 0{G Af) as expected, and 

Qoo = Rqo — 2^700-^ = ~3^^i' 
1 2 

Qij = Rij - -5ijR= --5ijGA\. (Ill) 

In the presence of some local matter distributions p(x), this background will be dominated 
by the gravitational field due to p(x). This is contained in R^^} which we compute next. 



Linear matter contribution. Now consider i?^^, which is the desired contribution lin- 
ear in p. To proceed we consider the time-and space-like components separately: 



i?oo = -2cu29%9a50 + 87rcu\2 



A\ 

2uj^goSg + Snu^glj^ ^ AnGp + 2U{x)uj^ 
Ai 



inGp (112) 
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where U{x) is the gravitational potential due to p as determined earlier; that term is 
negligible since AirGp ~ AU{x) S> U{x)u'^ by assumption. This is in agreement with GR 
as discussed below. Roi = follows again from do(j)di(j) = 0, and 

~ -{didkgo){djdkSg) - {djdkgo){didkSg) - Snuj'^j^godidjgQ 

^ '^u^god,d,5g+hnGp5,j (113) 



using (103). This gives 



R = -47rGp+'^uj^goASg + 4nGp = -'^4nGp (114) 



and 



1 2 

Qoo = Roo - 2^700^ = gSTrGp, 

1 2 2 2 

Qij = Rij - -SijR ^ -u godjdi6g + -iirGpSij. (115) 

In particular, in regions where p{x) = we obtain indeed 

Goo = Roo = GToo = 0, R = (116) 

as in GR, up to corrections due to the vacuum energy. Recall that in GR, the Einstein 
equations for = p, T*-' = imply Qqq = SnGp, and furthermore Rij = AngijGp, R 



00 



AnGp and SirGp = R. Thus (112) and (116) agree with GR, while the space-like (pressure) 
components of Rij and Qij are different here from GR. In particular, for p = the essential 
difference to GR is an anisotropy of Rij. Indeed we should not expect complete agreement 
with GR since there are no harmonic embeddings of non-trivial Ricci-flat 4-manifolds [37] . 

The non- vanishing components of Q"^^ = 0{Gp) would endanger the Newtonian limit if 
indeed A4 = Apianck- On the other hand, this does not pose any problem if A4 <^ Apianck, 
so that the contributions of A'^Q'^^ are negligible. This scenario is quite possible and in fact 
very appeahng as discussed below, greatly reducing the required fine-tuning for Ai. 



4.3 Quantization, scales and stability 

Perhaps the most remarkable result is that gravity arises in the harmonic branch of the 
matrix model, without even using the induced gravitational action i.e. the Einstein-Hilbert 



term. This means that the scale A4 in front of the induced Einstein-Hilbert term (34) could 
be much smaller than the Planck scale, 

A4 < Apianck (117) 
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because the Newton constant G is determined dynamically through (71). In principle, A4 



might be as low as 0{TeV) In fact this should be so in the present context, since 



otherwise the non- vanishing components Qij found above would enter equation (39) and 
might spoil the results of the previous sections (unless they vanish in a more sophisticated 
solution). This is of course a very attractive scenario, because then the difference of scales 
between Ai and A4 would be greatly reduced, largely resolving the cosmological constant 
problem. 

Another extremely interesting feature of the harmonic branch of the model is that the 
quantization of gravity seems to be straight-forward and well-behaved. Indeed A4 is essen- 
tially the scale of = 4 SUSY breaking, and the model can be viewed as non-commutative 
A^ = 4 SYM on which is expected to provide a well-defined quantum theory. The de- 
grees of freedom in 9^'^{x) can be interpreted as would-be U{1) non-commutative gauge 
fields on Mg, which are also well-behaved in the IKKT model. The excitations of the scalar 
fields (pi are essentially harmonic excitations due to the brane tension/ vacuum energy 

^vac = -2A.\ j rf^xVH ~ -Ai j d'x{2 + r]''''d^<j)'d,<j)'6,j) (118) 
with positive energy 

Evac ~ Af / d'x (2 + 9o05o0 + d,(l)di(l)6'^). (119) 



This is similar to the quantization of a Klein-Gordon field. Note in particular that the 
energy of the gravity bags is proportional to QqUj'^, hence its scale is more-or-less deter- 
mined dynamically by the initial energy distribution. In particular, it is plausible that the 
total energy of the universe consisting of brane fluctuations and the total mass is positive. 
This suggests that also after structure formation, the (vacuum) energy of the gravity bags 
associated to a galaxy with mass M should be no less than its gravitational binding energy 
MU. Then increasing would increase the total energy, which is important for stability 
reasons. 

In contrast, many of the difficulties associated with the quantization of GR are expected 
to be recovered in the Einstein branch of the model. 

It remains to be seen whether the emergent gravity is compatible with the precise solar 
system constraints, which requires a more complete analysis. Finally it should be pointed 
out that the considerations of this paper should be extended to cover the case of small 
extra dimensions, e.g. in the form of fuzzy sphere(s). The basic results of this paper are 
expected to apply also in that case. 

A similar possibility has also been considered in [9] however in a rather different context of higher- 
dimensional GR with branes. 
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4.4 Cosmological context and perturbations 



We briefly explain how the above solutions fit into a consistent cosmology. Assuming that 
the vacuum energy dominates the energy density due to matter, cosmological solutions 
of emergent NC gravity were obtained in [11] as harmonically embedded branes A4'^ C M}^ 



( 



\ 



I sinh(x) sin^cosy^ ^ \ 
sinh(x) sin^sin(y9 
sinh(x) cos 6' 
V cosh(x) 


Xc(t) 



(120) 



where 



Here ri^h = diag(+, ...+, - 
and has a FRW geometry 



7^(t) = a{t) 
+, +) for k = 



cos ip it) 

-1. This embedding is harmonic, AgX 
dT? = dx^ + sinh2(x)rf^]^ 



121) 
= 0, 

:i22) 



corresponding to spatial curvature k = —1. The physics of the early universe in this model 
is quite different from standard cosmology and requires a more detailed analysi^^ however, 
it is a solid prediction that the solution approaches a{t) — >■ t for late times, i.e. a Milne 
universe. Remarkably, this geometry is in good agreement with the basic observational 
data including the tye la supernovae data [38], which are usually interpreted in terms of 
an accelerating universe. 



Milne universe. The Milne geometry is expected to hold in this model as long as the 
vacuum energy A^ dominates the energy density due to matter. In particular, there is no 
upper bound on Ai from cosmology. These statements are easy to understand, recalling 
that the Milne universe is nothing but (a sector of) fiat Minkowski space M'^ C M^°. Indeed 
consider the fiat metric ds"^ = —dr"^ + dr"^ + r'^dD? on M^. In terms of the new variables 

r = tcosh(x), r = tsinh(x), (123) 

this metric takes the form of a FRW metric with a{t) = t and = — 1, 

ds] = -dr"" + dr^ + r^dn^ = ^ (-c/t^ + t^dx^ + smh\x)dn^)) 

= -dt^ + a{t)\dx^ + smh\x)d^^)- (124) 

In particular the consistency with the CMB data cannot be reliably addressed at this point, see 
however [38] for the simplified case of an exact Milne universe. The above refined solutions yield a big 
bounce and an early phase with power-law acceleration [11], but a proper treatment of matter is still 
missing. 
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Clearly, fiat Minkowski space is a solution of a brane with tension resp. vacuum energy 
~ Af, and in fact it is stabilized by a large Ai. 

We claim that matter such as stars and galaxies lead to local perturbations resp. fluctu- 
ations of this cosmological solution, and the previous results including Newtonian gravity 
and long-distance deviations go through with very minor adaptation. This can be seen 
easily by using the previous localized solutions on and rewriting them in terms of the 
Milne variables. For example, the spherically symmetric solutions centered at the origin 
r = X = become 

0o(x) = ^7(r)e- = ^7o^^e- 

our 

sinjutsinhix)) ^i^coshMt 
out smh(xj 

« J„?M^e"' (125) 

uja{t)x 

using X <C 1 in the neighborhood of the mass. Thus in comoving Milne coordinates, this 
looks like a spherically symmetric solution with red-shifted wavelength a{t)uj = too. Note 
that the effective gravitational constant G = °4 is independent of t, which is reassuring 
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A particularly interesting aspect is the slow "extrinsic" rotation (121) in the embedding 
Ai C M'^ for small t, which turns out to be [11] 

~ ^- (126) 

This might be very important in the context of fluctuation spectra in the early universe. 
Clearly a more complete perturbation analysis is required, in particular in order to address 
the issues of fluctuation spectra in the context of the cosmic microwave background. 

We can thus summarize the physical aspects of the model as follows, leaving aside its 
theoretical appeal e.g. with respect to quantization. The most attractive feature of the 
model is that it naturally predicts a (nearly-) flat universe, resulting in luminosity curves 
e.g. for type la supernovae which are close to the observed ones (usually interpreted in 
terms of cosmic acceleration) without any flne-tuning. On the other hand, reconciling 
it with observed gravity on scales less than or equal to galactic scales does impose an 
upper bound on the vacuum energy of order eV. This in turn might offer a mechanism 
for (partially?) explaining the galactic rotation curves, without requiring large amounts of 
dark matter. 



Gravitational (non—) constant G. Up to now, the parameters uj,go and therefore G 
were undetermined. However these are dynamical quantities, in particular they depend on 
the initial conditions. Thus we have to face the question why in particular G should be 

^^G might change in the very early universe due to deviations from the Milne metric. 
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the same in different parts of the universe, in particular in different galaxies. Indeed one 
should expect in this model that G depends somewhat on the individual galaxies, however 
the variation should certainly not be too large in our universe. 

While we cannot offer a completely satisfactory answer yet, some insight can be obtained 
from the above cosmological solutions. We have seen that deformations g{x) of the brane 
due to e.g. galaxies can be considered as perturbations propagating in an approximately flat 
Milne resp. Minkowski background, with uj and Qq remaining essentially constant in time. 
In reality, we know that galaxies are typically parts of larger structures such as (super)- 
clusters and filaments. This large-scale structure is expected to provide the dominant 
contribution to g{x) and therefore to G, which is therefore related over cosmological scales. 
Moreover, it is very plausible that the dynamics of structure formation in the universe 
leads to similar scales for gQ{x) of these dominant cosmic structures, given the high degree 
of homogeneity of the initial conditions seen in the CMB background. Furthermore, the 



rotation ip of the above cosmological solution (121) might provide a natural seed for the 
required rotation of the large-scale gravity bags. On the other hand, there could be some 
other more rigid mechanism for stabilizing G which we did not identify here. 



5 Discussion 



We obtained in this paper some remarkable results on the long-distance properties of emer- 
gent gravity in Yang-Mills matrix models, notably the IKKT model. Space-time is modeled 
by a 3+1-dimensional noncommutative brane solution, which acquires an "emergent" met- 
ric. Its dynamics is governed by the brane tension as well as additional gravitational terms 
including the Einstein- Hilbert term in the (quantum) effective action. There are two types 
of gravitational solution: an "Einstein branch" which is very similar to general relativity, 
and a "harmonic branch" where the branes are governed by a brane tension. We focus on 
the harmonic branch in this paper, and study the deformation of the space-time brane and 
its effective metric due to static localized mass distributions, having in mind e.g. galaxies. 

Due to the brane tension, the basic excitations of the space-time brane are essentially 
harmonic waves. Large matter clusters such as galaxies are embedded in such "gravity 
bags", which are rotating standing waves of the embedding with long wavelength L^. Stan- 
dard Newtonian gravity is recovered inside these "gravity bags" due to local matter such 
as stars. Moreover, there is an effective gravitational constant Aeg inside the bags, which 
leads to a significant enhancement of orbital velocities at large distances, quite reminiscent 
of the observed galactic rotation curves. At very large distances > L^, the harmonic em- 
bedding leads to a screening ~ ^ of the gravitational potential, reconciling relatively large 
vacuum energies (of order eV , say) with a consistent cosmology similar to a Milne universe. 
The latter is known to be in remarkably good agreement with the basic cosmological con- 
straints, leaving aside the CMB fluctuations which require a detailed understanding of the 
early universe in this model. 
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These results have important physical implications. An obvious conjecture is that the 
observed enhancement of the galactic rotation velocities compared with the Newtonian 
law is primarily due to the above results, i.e. an effective vacuum energy Aeff inside the 
gravity bags. This gets additional support by the result that the gravitational constant G 
and therefore Acfj is not universal but determined dynamically, and may therefore differ 
somewhat from galaxy to galaxy. This does not imply that there is no dark matter at 
all, but the model clearly requires far less dark matter in the galactic halos than what is 
invoked in the ACDM model. 

It is certainly remarkable how close this simple and rigid model comes to observation. 
In particular, large-scale cosmological observations seem to be reproduced much more nat- 
urally in the harmonic branch than in GR, without particular fine-tuning of the vacuum 
energy. The basic solar system observations can be reproduced to a good approximation, 
which seems to require an upper bound on Aeg on the order of eV . However, the preci- 
sion tests provide a challenge, and it remains to be seen whether they can be met in the 
harmonic branch of solutions under consideration. The only obvious modifications of the 
matrix model would be quadratic or cubic (soft) terms, which do not spoil the good UV 
behavior of its quantization. The cubic terms lead to compactified extra dimensions as 
fuzzy spheres, which is very natural and desirable for particle physics [39]. The quadratic 
terms might help to stabilizes the "scale" of space-time resp. the effective vacuum energy. 
Indeed observational constraints require that the vacuum energy is quite small (perhaps 
0{eV)). This may be due to a cancellation between the bare action and the induced 
quantum-mechanical vacuum energy, but quadratic terms in the matrix model might also 
play an important role here. 

Another remarkable feature of the model is that gravity is naturally weak. Indeed the 
scale of gravity is found to be G ~ ,4° 4 , where is a cosmological length scale, and Ai 

uj 1 

the vacuum energy. Moreover G can differ to some extent in different parts of the universe, 
and depends on the surrounding mass distribution which determines a "gravity bag" . 

In any case, an essential result of this paper is that Newtonian gravity arises in the 
matrix model simply due to the brane tension of a harmonically embedded space-time 
brane; the Einstein-Hilbert action is not required. This suggest the following very appeahng 
scenario: the coefficient A| of the induced Einstein-Hilbert term could be much sma//er than 
the Planck scale. This would also greatly alleviate the fine-tuning issue for the vacuum 
energy, and the precise form of the (induced) gravitational action is not essential. This 
basic mechanism might apply in a more general framework than that of Yang-Mills matrix 
models. However, the matrix model offers a clear concept of quantization, in particular 
the IKKT model. It contains not only (emergent) gravity but all ingredients required for 
a theory of fundamental interactions, in particular nonabelian gauge fields and fermions. 
This would provide a very interesting and accessible quantum theory of gravity, which has 
the potential to resolve some fundamental problems in this context. 

There are many obvious shortcomings in this paper which require more detailed work. 
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The present understanding is still at a rather crude level, and more detailed work is needed 
before claiming to seriously challenge the ACDM model. Probably the main problem are 
the space-like components of Q'^'^ resp. grr, which do not quite agree with GR even in 
vacuum. It remains to be seen if the solutions presented here can meet the solar system 
precision tests, or if a more complete solution taking into account notably 9^'^{x) will lead 
to the suitable corrections. Furthermore, the short-distance properties of a point-mass 
solution i.e. black holes should be studied in the matrix model framework. Here higher- 
order corrections due to 9'^^ are expected to play an important role. Finally, the Einstein 
branch of the IKKT model should provide a quantization of more conventional general 
relativity coupled to matter. 
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